9:371-386, 1970), X n is given by: 



16C 



(12) 



where Rei is the Reynolds number of the parent branch, C is a constant, D x and D n are the 
diameters for the parent branch and nth daughter branch respectively, and L n is the entrance 
5 length for the nth daughter branch, 
pi Values of the kinetic energy and dissipation factors are dependent on the nature of flow at 

^ the bifurcation. These loss coefficients are approximated by Wolf (Wolf, T. "An 
« Experimental/Theoretical Investigation of Parallel Inhomogeneities in Respiratory Flows." 
i J Ph.D. thesis, Dept. of Mechanical Engineering, Massachusetts Institute of Technology: June, 
40 1990; incorporated herein by reference) as: 



( i A 2 ( 1 ~\ 2 



1- 1 



V K c,n J 



and f n = 



_1 

\ K c,n J 



(13), (14) 



The contraction coefficient K c n represents the ratio between the minimum normal cross- 
sectional area of the streamtube within the separated region of the nth. daughter branch, and the 

15 area of the branch itself. Hence, for smaller angles between adjacent daughter branches, 
separation does not occur and K cn should approach unity. Using literature values for the 
anatomical branching angles of the arterial system, a linear relationship between angle and K c n 
is assumed and values assigned for each branch. The method is imprecise, however a sensitivity 
analysis of the contraction coefficients demonstrates that exact values are not required due to 

20 their small effect on the system as a whole. The energy conservation equation [Eq. (10)], and the 
unsteady form of continuity [Eq. (11)] are solved by the MacCormack predictor-corrector 
computational scheme(Amsterdam, E.A,, Banas, J., Cartley, J.M., et al. Clinical assessment of 
external pressure circulatory assistance in acute myocardial infarction. Am. J. Cardiol, 45:349, 
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1990; incorporated herein by reference), coupling the state variables of the bifurcation end-nodes 
with the internal points of the involved arterial elements. 

Ventricle . A model of the left ventricle acts as an upstream boundary for the arterial tree. The 
5 ventricle may be approximated as a chamber whose compliance (or the inverse of elastance) 

changes with time, thus driving flow through a unidirectional exit valve into the aorta. 

Following the work of Suga and Sagawa (Suga, H. and Sagawa, K. Instantaneous pressure- 
5 volume relationships and their ratio in the excised, supported canine left ventricle. Circulation 
M= Research, 35, 1974; incorporated herein by reference), the specified elastance curve E(t) is used 
■;1;0 to characterize ventricular ejection and filling. Extensive studies of the pressure-volume 

relationship in canine ventricles (Suga, H. and Sagawa, K. Instantaneous pressure- volume 
U\ relationships and their ratio in the excised, supported canine left ventricle. Circulation Research, 
y 35, 1974; incorporated herein by reference) have shown that the basic shape of the systolic 

portion of the pressure-volume curve remains unchanged, regardless of loading or wall 
15 compliance changes. Thus, the systolic wall elastance may be characterized by only two 

parameters; the peak wall elastance, E mQX , and time to this peak elastance during one cycle, 

T m8LX . The duration of the cycle is extrapolated from Suga and Sagawa' s data suggesting that 

?max spans approximately 30 - 50% of the total cycle. 

In a later extension of this model by Suga and Sagawa (Suga, H., Sagawa, K., Demer, L. 
20 Determinants of Instantaneous Pressure in Canine Left Ventricle: Time and Volume 

Specification. Circ. Res. 46: 256-263, 1980; incorporated herein by reference), the effects of 

myocardial viscoelasticity were taken into account by the addition of one term to the ventricular 

pressure-volume relation. We found this approach to better fit data for humans (Ozawa, E.T. "A 
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numerical model of the cardiovascular system for clinical assessment of the hemodynamic state." 
Thesis, Ph.D., Massachusetts Institute of Technology, September, 1996; incorporated herein by 
reference) and have consequently used the following form in the present simulations: 



The coefficient a is a scaling factor for the time-dependent viscoelastic effects, P vent is the 
ventricular volume, V vent is ventricular volume, V venUo is the zero-pressure filling volume, 
P is the transpulmonary pressure in the chest cavity which is expected to alter the left 



assume it to be a half sinusoid, whose duration and amplitude can be modified to represent 
j different compliances and heart rates. Thus, the ventricular pressure may be solved for, given 
j that the flowrate at the root of the aorta is equal to the time derivative of ventricular volume. 

15 Sinuses of Valsalva . As flow begins to reverse, the valve leaflets are swept backwards and close 
without sustaining a significant pressure gradient. Filling of the sinuses continues until the valve 
leaflets are maximally distended, at which time the leaflets are able to sustain a pressure 
gradient. This may be modeled as an abrupt decrease in the aortic root compliance to a new 
value, C sinus , as well as imposing a zero-flow boundary condition at the first node. This 

20 condition is held until the ventricular pressure again exceeds aortic pressure. 



5 




ia?0 ventricular transmural pressure 



:. We define the isovolumetric contraction curve £*(/) and 



Terminal branch points . The numerical model described here uses linear segments to represent 
the larger vessels in the main arterial tree, but modeling the finer branching structure 
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approaching the arterioles in this manner is impractical. Rather, the terminal vessels are modeled 
as a lumped parameter Windkessel (Berger ? D.S., Li, J.K-J, and Noordergraf, A. Arterial Wave 
Propagation Phenomena, Ventricular Work, and Power Dissipation. Ann. BiomeJ Eng. 23:804- 
81 1, 1995; Berger, D.S., Li, J.K-J, and Noordergraf, A. Differential effects of wave reflections 
5 and peripheral resistance on aortic blood pressure: a model-based study. Am. J. Physiol 266: 
(Heart Circ. Physiol) 35:H1626-H1642, 1994; each of which is incorporated herein by 
reference). The model allows the behavior of the small arterial vessel beds to be captured using 

jfi only a few parameters and accurately mimics peripheral wave reflections (Ozawa, E.T. "A 

numerical model of the cardiovascular system for clinical assessment of the hemodynamic state." 

rife Thesis, Ph.D., Massachusetts Institute of Technology, September, 1996; incorporated herein by 
reference). The Windkessel consists of a resistance R s in parallel with a compliance C s , where 
the resistance represents the pressure drop associated with the terminal arterioles, and the 
compliance represents the total compliance of the small artery network. In series upstream from 
the Windkessel is an additional element Z 0 , which represents the entrance impedance of the 
15 small arterial bed with an associated pressure drop of P-P c , where P is the pressure at the last 
node in the terminating arterial segment, and P c is the capillary bed pressure. This impedance is 
matched to that of the adjoining element to avoid the generation of reflections resulting from an 
impedance mismatch at this interface. Thus, Z 0 is approximated as pc I A using the wavespeed 
and area for the element node bordering the Windkessel. From the electrical analogue, the 

20 following equations may be written as a function of the resistances and capacitances: 



P -P dP 

_C * + C _JttL = £ (16) 

K s dt 
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Q = ±(P-P e ) d7) 



where Q is flow entering the Windkessel, P v is the venous pressure, and P tm =P C -P e is the 
transmural pressure, where the external pressure P e may be specified. The Windkessel is 
coupled numerically to its upstream element in this model using the method of characteristics. 
5 Within the region of external compression the externally applied pressure is used as the 

reference pressure for the capacitor in the windkessel model. Venous pressure is assumed 

5 constant for the purpose of these calculations. Consequently, the dynamics occurring on the 
£ venous side of the circulation, while potentially important, are not considered in the present 
W calculations. The implications of this assumption are discussed later. 

iio The model, as implemented for these simulations, is operated in "open mode" in that the 

6 venous and pulmonary circulations have been omitted. In so doing, the dynamics of the venous 
- "H bed associated with EECP are essentially ignored on the assumption that the changes in mean 
R venous pressure due to EECP will have minimal effect on the pulsatile flows and pressures on 

the arterial side. Variations of venous pressure in the region of compression may have a 
15 somewhat greater effect as discussed below. 

Parameter Specification . Information on cardiovascular system parameters exists, but the data 
are highly variable and extremely dependent on the state (e.g. posture) of the individual or 
animal at the time when measurements are taken. Nonetheless, a standard case that produces a 
20 reasonable model output for a given set of parameters is required. Avolio (Avolio, A.P. Multi- 
branched model of the human arterial system. Med. Biol. Eng. Comput. 18:709-18, 1980; 
incorporated herein by reference) presents an extensive list of arterial lengths, diameters, wall 
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thicknesses, and Young's moduli for most of the major human arteries (a sum total of 128). For 
the present model, 30 elements are used. This network is shown schematically in Fig. 1. The 
numbered elements correspond to major arteries whose properties are provided in Table 2. 



Table 2. 

sectional 
vessel. E 



Specifications for the 30 element model: arterial properties. Proximal and distal cross- 
areas are given, and it is assumed that area varies linearly with distance along the 
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1 internal carotid 
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20 


1 external carotid 


0.113 


0.196 


0.071 


8.0 


21 


1 radial 


0.229 


0.08 


0.08 


8.0 


22 


1 ulnar 


0.232 


0.139 


0.113 


8.0 


23 


coeliac 
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0.478 


0.478 


4.0 


24 


renal 
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0.212 


0.212 


4.0 


25 


sup mesenteric 


0.054 


0.581 


0.581 


4.0 


26 


inf mesenteric 


0.045 


0.08 


0.08 


4.0 


27 


profundis 


0.121 


0.166 


0.166 


16.0 


28 


post tibial 


0.306 


0.102 


0.102 


16.0 


29 


ant tibial 


0.295 


0.031 


0.031 


16.0 


30 


peroneal 


0.313 


0.053 


0.053 


16.0 
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Arterial elements were specified in the model by a proximal and distal internal radius, 
from which the cross-sectional areas were calculated. The area was assumed to be a linear 
function of length between bifurcations. All elements were discretized into nodes, separated by a 
5 spatial increment, nominally 0.01 m. The branching pattern of the arteries was also taken from 
the arterial tree layout given by Avolio (Avolio, A.P. Multi-branched model of the human 
arterial system. Med. Biol Eng. Comput 18:709-18, 1980; incorporated herein by reference). In 
i|l addition to geometrical measurements, the elasticity of each artery segment is specified in order 
K* to calculate the nominal reference wavespeed Co for each element, using the Moens-Korteweg 
;f|0 equation: 



I Here, c 0 depends locally upon Young's modulus E, inner radius R, fluid density p and wall 
15 thickness h. Both R and h were assumed known with the internal static reference pressure P 0 
equal to 100 mm Hg (13.3 kPa). Assuming a linear relationship between wall thickness and 
vessel radius, the model calculates c 0 at each node given dimensions and material properties. 
Values for Young's modulus (Table 2) were also obtained from Avolio' s original data. 



20 complicating the task of defining a "standard" distribution of flow. An alternative approach used 
in this model is based on the observation that the caliber of a vessel is related to the flow rate it 
conveys. This was first proposed by Murray in 1926 and later shown to result from biological 
factors that tend to maintain a constant wall shear stress in all arteries (Kamiya, A., Bukhari, R., 




(18) 



Local regulation of flow is extremely transient and dependent on a variety of factors, 
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and Togawa, T. Adaptive regulation of wall shear stress optimizing vascular tree function. Bull 
Math Biol 46: 127-137, 1984; incorporated herein by reference). Thus, assuming a blood 
viscosity of 4.0xl0~ 5 kg m" 1 s"\ and a mean wall shear stress of 1.5 N in 2 at a mean arterial 
pressure of 13.3 kPa (100 mm Hg), the required flow can be calculated for each node within the 
network. The difference in flow between nodes is assumed due Eq. (4). From the calculated 
required flow at each distal boundary node, the appropriate value of the terminal Windkessel 
resistance R s is determined for each terminal element. This method provides values of mean 
flowrate that are in reasonable agreement with known physiological values. Where the predicted 
values differ from literature values, adjustments have been made. Values of other hemodynamic 
parameters are shown in Table 3. 



Table 3. External Pressurization Input Control Parameters. All times referenced to the 
beginning of cardiac systole. 



Parameter 


Description 


Baseline Values 


Tinfl 


Time from onset of cardiac cycle at which pressure starts 
to rise 


0.20 s 


Tdefl 


Time to begin deflation of cuffs from maximum external 
pressure 


0.72 s 


T ca rd 


Period of cardiac cycle 


0.86 s 


Pcalf 


Maximum external pressure applied to calf vessels 


200 mmHg 


Pth 


Maximum external pressure applied to thigh vessels 


150 mmHg 


Pla 


Maximum external pressure applied to lower abdomen 
vessels 


100 mmHg 


At seg 


Time interval between inflations of adjacent cuff regions 
with the proximal region always pressurized first 


0.03 s 


tramp 


Time it takes pressure to rise to and fall from its maximum 
value 


0.03 s 


APseg 


Pressure difference between cuffs (pressure always 
increasing in the direction of the foot) 


50 mmHg 


P m 


Average applied pressure in the three cuff regions 


150 mmHg 
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A base state was chosen, typical of conditions used clinically, from which the effects of 
various parameter variations could be studied. All data presented are taken from the tenth heart 
cycle of the model to ensure that the simulation has reached a steady state. A heart rate of 72 
5 beats/min is used for all simulations. Values of the control parameters used for this base state are 
given in Table 3. Some judgment was exercised in parameter selection. A second set of 
parameter values, with peak ventricular contractility reduced from 6000 to 1000 dyn/cm 5 , 
yQ systemic vascular resistance increased from 1000 to 2666 dyn/cm 5 and end diastolic volume 
\f increased from 120 to 280 ml was used to simulate a patient with compromised ventricular 
;;j0 function. Validation of the model with normal parameters included comparisons to measured 
J" waveforms at various locations in the arterial system and measured arterial input impedance, 
hi Comparisons were also made to pressure and velocity traces found in the literature, but since 
bj these vary considerably among individuals direct comparisons are of limited value. These can be 
H ! found in Ozawa (Ozawa, E.T. "A numerical model of the cardiovascular system for clinical 
15 assessment of the hemodynamic state." Thesis, Ph.D., Massachusetts Institute of Technology, 
September, 1996; incorporated herein by reference). 

External Pressurization Scheme. A three-step graded-sequential compression procedure was 
employed in all the simulations presented here. In sequential compression, a wave of 
20 compression is applied to the vessels by inflating the three pressurization cuffs for the calves, 
thighs, and lower abdomen sequentially from ankle to groin. The pressure level applied by the 
cuffs decreases from calf to thigh, and from thigh to lower abdomen cuffs. In contrast to the 
emptying behavior characteristic of uniform compression, sequential compression produces a 
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collapse in the vessels that proceeds from the foot toward the heart. Thus, the blood is 
effectively "milked" from the vessels in the lower extremities and does not pass through a 
constrictive throat as in uniform compression (Lueptow, R.M., Karlen, J.M., Kamm, R.D., 
Shapiro, A.H. Circulatory Model Studies of External Cardiac Assist by Counterpulsation. 
5 Cardiovascular Research, 15(8):443-455, 1981; incorporated herein by reference). In graded 
compression the maximum level of pressure attained in each segment is greatest in the periphery 
and falls in the direction of the heart. The application of graded compression also helps to 
C eliminate the occlusive throat and, in combination with sequential pressure application, produces 
fZ rapid and complete emptying of the vessels (Lueptow, R.M., Karlen, J.M., Kamm, R.D., Shapiro, 
5t0 A.H. Circulatory Model Studies of External Cardiac Assist by Counterpulsation. 
7 Cardiovascular Research, 15(8):443-455, 1981; Zheng, Z.S., Li, T.M., Kambic H., et al 
yj Sequential external counterpulsation (SECP) in China. Transactions of the American Society of 
W Artificial Internal Organs, 29:599-603, 1983; each of which is incorporated herein by reference). 
^ The cuffs used to provide pressurization of the lower extremities in EECP are modeled as 

15 external pressure sources on the lower abdomen, thigh, and calf arteries. To simulate graded- 
sequential compression in the model, the arterial tree elements for the lower body are divided 
into three regions, shown in Fig. 1, representing the areas covered by the three pressurization 
cuffs in EECP. 



20 External Pressurization Control Parameters. Clinical and computational studies have shown the 
efficacy of EECP depends upon the mode of operation and parameter values used to control the 
device (Bai, L, Wu, D., Zhang, J. A Simulation Study of External Counterpulsation. Comput 
Biol Med.,. 24(2): 145-156, 1994; incorporated herein by reference). These parameters include 

3184967_1.DOC 36 of 60 



the cuff inflation and deflation timings, the maximum pressure level applied externally to the 
vessels by each cuff, and the time delay of pressurization and depressurization between the calf, 
thigh, and lower abdomen cuffs for sequential compression. Table 3 shows a detailed 
description of the individual input control parameters governing external pressurization in the 
5 model. 

In clinical practice, the application of external pressure during EECP is timed with the 
patient's electrocardiogram. In the EECP model, this process is accomplished by adjusting the 
I: timing of applied external pressure in each of the three compartments relative to left ventricular 
I- contraction, as characterized by E(t). For graded-sequential compression, the pressure in each 
JffO cuff rises linearly to its maximum value over a time t^, is held constant until a time T defl , and 
; M then falls linearly over a time t^. The calf, thigh, and lower abdomen cuffs are inflated at times 
R T. T. +At andT +2At respectively. The maximum applied pressure is decreased between 

1 7. ?. infl' infl seg' mfl seg xr 

I'?: the calf and thigh cuffs and the thigh and lower abdomen cuffs as specified by P calf , P.,, and P la . 
lI The cuff deflation time, T def) , is the same for all three cuffs to simplify the parameter study. The 
1 5 parameters used in the temporal application of external pressure during the heart cycle are given 
in Table 3. For all other parameter values, see Ozawa (Ozawa, E.T. "A numerical model of the 
cardiovascular system for clinical assessment of the hemodynamic state." Thesis, Ph.D., 
Massachusetts Institute of Technology, September, 1996; incorporated herein by reference). 

Mean applied pressure was chosen at a level thought to produce minimum trauma to the 
20 patient while still providing a reasonable measure of benefit. The pressure increment between 
segments was viewed as sufficient to prevent proximal arterial collapse and a consequent 
impairment of vessel emptying while still providing ample pressure at the lower abdomen region 
to produce significant emptying. Consistent with the notion that arterial emptying should 
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proceed at a speed comparable to the speed of the arterial pressure pulse (about 8 m/s in the 
peripheral arteries), the time delay between segment compressions was chosen to be 
approximately equal to the wave transit time through each of the pressurized compartments. 
Pressure rise time, as shown by Bai et aL, should be as short as possible (Bai, J., Ying, K., Jaron, 
5 D. Cardiovascular responses to external counterpulsation: a computer simulation. Med. Biol 
Eng. Compute 30:317-323, 1992; incorporated herein by reference). Therefore, a value was 
chosen close to the practical lower limit. 

External Pressurization Measures of Merit. The effectiveness of enhanced external 
y l0 counterpulsation is assessed in terms of the following measures of merit: 

J;i Mean Diastolic Pressure. The increase in diastolic pressure, or diastolic augmentation, is 
\ J] characterized by the mean diastolic pressure ratio: 




where P a0 rtic is pressure at the aortic root, T D and T s are the times at which diastole and systole 
end, and the subscripts "compr" and "0" refer to cases with and without external compression, 
respectively MDP is an indication of how diastolic pressure is increased with pressurization. All 
20 pressures are measured when the model has reached steady-state after 10 heart cycles. 
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Mean Systolic Pressure. The effect of EECP on systolic pressure is quantified using the mean 
systolic pressure ratio: 



MSP = 



rjj Jq aortic 

Js 


0 


-Ifp .dt 

T Jo aortic 
_ 1 S 


compr 




±\ Ts P dt 

rp Jq aortic 

_ s 


0 



MSP is a measure of the extent of left ventricular afterload reduction with pressurization. 



Emptying Effectiveness. The blood volume emptied from the leg vessels enters the aorta, and 
hence determines the amount of diastolic augmentation achieved by EECP. It also provides a 
measure of the extent to which arterial diameter changes, and therefore relates to vessel wall 
strain. The emptying effectiveness parameter, EE, is used to measure the efficiency of the 
emptying process for the vessels receiving pressurization. EE is calculated for a single vessel 
using the equation: 



\\a dx\-\\A dx] 

EE= L -, *" mpr (21) 



where A is the cross-sectional area of the artery. The integrations are taken over the entire 
region of pressurization for the artery of interest. For the "compr" case, arterial area is measured 
at maximum pressurization in diastole just prior to cuff deflation. The arterial area for case "0" 
is measured at the time step just preceding pressurization. Thus, the emptying effectiveness of 
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the artery represents the extent of arterial collapse under maximum pressurization with respect to 
the state of the artery just prior to pressurization. The state of the artery prior to pressurization is 
considered since the artery may be partially collapsed if there has not been sufficient time for it 
to completely refill. 

Shear Stress Index. An approximate measure of shear stress is defined, that accounts for the 
changes in cross-sectional area and flow velocity that accompany EECP. In the case of steady, 
fully-developed, laminar flow through a vessel of circular cross-section, wall shear stress could 
be computed as follows: 



(22) 



' where V is the mean flow velocity. Recognizing that as an artery collapses, its cross-section 

will likely deviate from circular, and that the flow is clearly not fully-developed nor steady, we 
1 5 will still assume to a rough approximation that 



T « «— )= (23) 



Actual values of shear stress will be larger than this due to the change in vessel shape and 
20 unsteadiness. However, as these effects are difficult to estimate accurately without resorting to 
fully three-dimensional calculations, we have chosen instead to use eq. (23) for the purpose of 
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estimating the relative values of shear stress between simulations with and without EECP. 
Accordingly, a shear index, S, is defined as 

3 3 
/ j w, compr / a *w,0 

S = ^ 3 ^ (24) 

5 

S where r w compr is the time integral of the wall shear stress for one cycle evaluated at the mid-point 

of the compression zone and T w 0 is the same value with no external pressurization. The 
% I summation sign indicates that the values of t w compr and t w 0 are summed over the three 
g compression zones. 

uo 

W Results 

f " Pressure pulses at the radial artery and aortic root computed by the model (Fig. 3) with 

and without graded-sequential external compression from the lower abdomen to the foot clearly 
illustrate the hemodynamic effects of EECP. Compression of a "normal" subject (Fig. 3a) is 
15 contrasted to EECP in a patient with reduced ventricular function (Fig. 3b). In both cases, 
pressure is applied by a three-compartment cuff with maximum pressures of 200, 150 and 100 
mmHg along the lower leg, upper leg, and lower abdomen, respectively. Note that in this 
instance of arterial counterpulsation (pressure application during cardiac diastole and release of 
pressure during systole) systolic pressure is reduced while diastolic pressure is augmented, 
20 leading to the combined effects of reduced ventricular afterload and enhanced coronary blood 
flow. 
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The effectiveness of EECP can be viewed in terms of the measures of merit defined 
previously. These results, computed for the "normal" subject, are shown in Table 4. Numbers 
shown in the table correspond to the fractional change in each measure from the case without 
compression and are defined so that they range in value from zero (no effect) to order one. 

5 

Table 4. Values for each of the measures of merit for the baseline conditions given in Table 3. 





Values under baseline conditions 


MDP 


0.0782 


MSP 


0.0238 


EE 


0373 


S 


3.22 



fll For this same condition of external compression, the time-varying arterial cross-sectional 

f_ area and a measure proportional to the time-varying shear stress (see Eq. (23)) are plotted for 
•JiO three locations (lower abdomen, thigh, and calf) in Figs. 4 and 5. During pressure application, 
[ the arteries collapse with sufficient speed to cause a flow reversal throughout much of the arterial 
network and a significant increase in vascular shear stress in the arteries of the lower extremity. 
The arteries in the lower abdomen and thigh (Figs. 4(a) and 4(b), respectively) refill rapidly upon 
pressure release, even rising to slightly above normal levels due to the strong compression wave 
15 generated and its reflection from the peripheral vascular bed. During refilling, shear stress 

attains levels roughly 3- to 4-fold higher than under normal conditions at all three locations (Fig. 
5). Features of particular interest in the context of endothelial function are the high shear 
stresses of reversing sign, and the significant arterial wall strain due to arterial collapse. 

One way to elucidate the relative effects of the various compression parameters is by the 
20 sensitivity matrix of Table 5. 
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Table 5. Sensitivity matrix. Each numerical value represents the fractional change in the 
particular measure of merit [A(MM)/MM m&an ] divided by the fractional change in the parameter 
[AY/Y mean ] as defined in Eq. (25). Note that a positive value indicates an increase in magnitude 



Y = 


P m 


AP se g 


tramp 


Atgeg 


Tinfl 


MDP 


1.01 


-0.140 


-0.107 


-0.086 


-0.887 


MSP 


1.26 


-0.512 


-0.147 


0.067 


-0.059 


EE 


-0.614 


0.236 


-0.075 


-0.029 


0.702 


S 


2.82 


-0.193 


-0.028 


0.107 


0.616 



5 Considering that each measure of merit (MM) (e.g., mean diastolic pressure, MDP) is a function 

h of each of the adjustable parameters (Y) (e.g., mean applied pressure), then an entry in the table 

m (X) represents the change in the measure of merit [A(MM)] divided by the fractional change in 

:;30 the parameter value: 



\A(MM)/MM me J 



(25) 



Large values indicate strong correlations between the measure of merit and the particular 
15 parameter. For example, a 10% increase in mean applied pressure (P m ) will produce an increase 
of 0.0078 (=0.1 1.01 0.0782) in MDP. This table is useful for identifying the parameters that, 
when varied, will have the greatest influence on the measure of interest. 

While the dependencies are generally quite symmetric about the baseline case suggesting 
a nearly linear dependence, there is one notable exception. The time to initiate cuff inflation 
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(Tjnfi) selected for the baseline case was near the optimum in terms of reducing mean systolic 
pressure. Thus, although T infl has a strong influence on most of the measures of merit, increasing 
or decreasing it by small amounts from the baseline case does not appear to have much effect on 
MSP. Changes of larger magnitude, however, will have significant deliterious effects; e.g., if 
5 cuff inflation occurs earlier, before the end of systole, systolic unloading will be severely 
affected. 

Mean applied pressure clearly has the greatest potential to enhance diastolic pressure and 
h2 increase levels of shear stress, although increasing mean pressure probably has the largest 

negative impact on patient tolerance. Altering cuff inflation time also exerts an important 
Z jo influence, although some of its effects are counter-productive (e.g. , when S and EE increase, 
; lf MDP falls). Reducing pressure rise time (t ramp ) is also beneficial, although this may be difficult 

to accomplish in practice. 

1 ™ Discussion 

15 The simulations of EECP presented here, provide insight into the dynamic processes that 

accompany EECP. Beginning near the end of systole, compression produces arterial collapse, 
sending a wave of retrograde flow up into the aorta, increasing pressure up as far as the aortic 
root and, presumably, augmenting coronary blood flow. When compression is released near the 
end of diastole, the arteries begin to refill, initiating a rarefaction wave that propagates toward 

20 the heart, reaching it at a time that produces systolic unloading. The extent of emptying of the 
leg arteries decreases toward the periphery, but corresponds to approximately half their normal 
arterial volume within the range of pressures tested here. Refilling to normal volumes is 
achieved in the most proximal arteries, but even at the levels of pressure used in these 
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simulations, is incomplete in the lower leg. Increasing pressure applied to the calf up to 300 
mmHg (results not shown) compromises refilling even further. The results in Table 3 
demonstrate the potential for controlling mechanical events related to cardiac assist or vascular 
cell stimulation, and perhaps more importantly, show that optimization of the measures relating 
5 to cardiac function is not always consistent with optimization of the vascular stimulus. 

In previous models of EECP, the arteries were represented as a collection of lumped 
elements and were therefore not capable of accurately capturing many of the phenomena 
3 associated with wave propagation through the arterial network and arterial collapse (Bai, J., Wu, 
h D., Zhang, J. A Simulation Study of External Counterpulsation. Comput Biol Med.,. 
U jO 24(2): 145- 156, 1994; Bai, J. ? Ying, K., Jaron, D. Cardiovascular responses to external 
s ^ counterpulsation: a computer simulation. Med. Biol Eng. Comput, 30:317-323, 1992; each of 
H which is incorporated herein by reference). The present model, though still discretized, solves 
JH the distributed differential equations and also incorporates the nonlinearities associated with 
m- arterial collapse and convective acceleration which are critical under conditions of EECP. 
15 Consequently, the model captures the influence of forward and backward propagating waves, can 
reproduce the complex impedance of real arterial networks (Ozawa, E.T. "A numerical model of 
the cardiovascular system for clinical assessment of the hemodynamic state." Thesis, Ph.D., 
Massachusetts Institute of Technology, September, 1996; incorporated herein by reference), and 
thereby provides a means to examine the detailed flow dynamics associated with EECP. As seen 
20 in Fig. 4, the onset of compression at the lower leg sends a surge of blood toward the heart, 
producing the rapid rise in cross-sectional area in the thigh (see e.g., t=3.60 s in Fig. 4(b)) and 
the lower abdomen (t=3.64 s in Fig. 4(a)) just prior to compression of these regions. The abrupt 
fall in cross-sectional area in the calf (beginning at t=3.56 s, Fig. 4(c)) is followed by an equally 
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abrupt rise in area, before the area decreases more consistently as pressure is maintained. The 
low frequency oscillation occurring in the iliac artery (t=2.9-3.2 s, Fig. 4(a)) is evidence of wave 
reflection from the proximal end of the aorta, causing some refilling while pressure is still 
maintained. These waves are highly damped, however, and are not seen in the thigh or calf 
5 regions. 

These results, in terms of the magnitude of the effect observed in arterial blood flow and 
pressure, can to some extent be compared with previous observations. In the multi-center study, 
y3 the hemodynamic effects of EECP were monitored by determining an "effectiveness ratio" 
H defined as peak diastolic pressure minus end-diastolic pressure divided by peak systolic pressure 
;P0 minus end-diastolic pressure (Arora, R.R., Chou, T.M., Jain, D., Fleishman, B., Crawford, L., 
7* McKiernan, T., Nesto, R.W. The multicenter study of enhanced external counterpulsation 
HI (MUST-EECP): effect of EECP on exercise-induced myocardial ischemia and anginal episodes, 
y Am. Col. Cardiol, 33(7): 1833-1840, 1999; incorporated herein by reference). Using the 

H= aortic blood pressure trace, this ratio is 0.90 for our standard case (Fig. 3b) and 1.60 for the case 
15 with a compromised ventricle ("diseased") simulation (Fig. 3d) compared to an average of 1.41 ± 
0.51 in the multi-center trial. Note that the patient values (1) were based on measurements with 
a finger plethysmograph and therefore are not directly comparable to our predictions, and (2) 
were obtained using maximum pressures up to 300 mmHg. In a separate study, most effects 
(e.g., change in cardiac output, ratio of retrograde to antegrade aortic flow) had nearly reached 
20 their maximal effect when the diastolic-to-systolic pressure ratio reached values in the range of 
one to 2 (Suresh, K., Simandl, S., Lawson, W.E., Hui, J.C., Lillis, O., Burger, L., Guo, T., Cohn, 
P.F. Maximizing the hemodynamic benefit of enhanced external counterpulsation. Clinical 
Cardiology, 21(9):649-653,1998; incorporated herein by reference). 
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We chose to study a range of pressures below those currently used clinically in 
recognition of the relatively high number of adverse experiences reported by patients receiving 
EECP. In the multi-center study (Arora, R.R., Chou, T.M., Jain, D., Fleishman, B., Crawford, 
McKiernan, T. ? Nesto, R.W. The multicenter study of enhanced external counterpulsation 
5 (MUST-EECP): effect of EECP on exercise-induced myocardial ischemia and anginal episodes. 
Am. Col Cardiol, 33(7): 1833-1840, 1999; incorporated herein by reference), 54.9% of 
patients experienced adverse effects, with the majority of these being device-related. The 
J: number of device-related adverse effects was reduced nearly 4-fold (from 37 to 10) in a separate 
'{?•: group of patients in whom applied pressures were decreased from 300 mmHg to 75 mmHg. Our 
rid results indicate that significant hemodynamic effects, especially in terms of enhancing arterial 
i shear stress and arterial wall strain, can be achieved with the use of considerably lower pressures, 
^ with mean values in the range of 150 mmHg. 

"if Our assumption of constant venous pressure has several potential implications. Consider 

first the effect of compression on the veins of the lower extremity. The magnitude and frequency 

15 of compression will almost certainly cause these veins to be in a state of collapse throughout the 
cycle since venous valves prevent rapid refilling as occurs on the arterial side. Thus, during 
compression, venous pressure will be at or slightly elevated above central venous pressure, but 
during the relaxation phase it will be strongly negative. We anticipate that this would lead to an 
increase in mean limb flow but of unknown magnitude. A second effect of compression, 

20 secondary to the collapse of the veins of the lower extremities, would be a slight elevation of 
central venous pressure leading to enhanced venous return and activation of the atrial 
baroreceptors. The effects of these changes are difficult to predict with the present model, and 
deserve further attention. For the purpose of predicting the pulsatile changes in flow and cross- 

3184967_1.DOC 47 of g Q 



sectional area in the leg arteries, however, these changes on the venous side should have 
relatively little impact. Several additional simulations (results not shown) were conducted with 
end diastolic volume increased by 20 ml, and with venous pressure increased 5 mmHg. The 
latter produced results indistinguishable from those presented here. The increase in end diastolic 
5 volume led to an overall increase in arterial pressure by about 10 mmHg, but the pressure profile 
with EECP changed very little. The net effect would be a reduction in systolic unloading and an 
increase in diastolic augmentation. 
^ Another potential source of uncertainty relates to the values for elastic modulus used for 

!! the arterial tree. Although we used the only data we were aware of in the literature 3 , these are 
; j0 only estimates and are therefore subject to error. Since the stiffness of the arteries under collapse 

is based on these estimates, it is possible that the pressures required to produce a certain degree 
jj4 of arterial emptying may also be uncertain. Add to this the subject-to-subject variability likely to 
l AJ be present, and it is apparent that these results should be used only as a rough guide to estimating 
; the parameter values for any particular subject, and that some amount of empirical testing is 
1 5 essential in practice. 

It is interesting to note that recent clinical results show benefits in cardiac function from 
as little as one hour of treatment per day. Under this protocol, Lawson et al, found that 17 out of 
18 patients receiving EECP for as little as 36 one-hour treatments reported improvement in 
anginal symptoms, despite prior medical and surgical therapy (Lawson, W.E., Hui, J.C, Soroff 
20 H.S., Zheng, Z., et al. Efficacy of enhanced external counterpulsation in the treatment of angina 
pectoris. American Journal of Cardiology, 70(9):859-862, 1992; incorporated herein by 
reference). The results of the recent multi-center study involving 139 patients confirmed these 
findings and showed that EECP reduces exercise-induced ischemia and angina (Arora, R.R., 

3184967_1.DOC 48 of 60 



Chou, T.M., Jain, D., Fleishman, B., Crawford, L., McKiernan, T., Nesto, R.W. The multicenter 
study of enhanced external counterpulsation (MUST-EECP): effect of EECP on exercise-induced 
myocardial ischemia and anginal episodes. J. Am. Col Cardiol, 33(7): 1833- 1840, 1999; 
incorporated herein by reference). A 1996 study by Lawson et al. also showed improved 
exercise tolerance in 22 out of 27 patients with chronic stable angina (Lawson, W.E., Hui, J.C., 
Zheng, Z.S., Burgen, L., Jiang, L., Lillis, O., Oster, Z., Soroff, FL, Cohn, P. Improved exercise 
tolerance following enhanced external counterpulsation: cardiac or peripheral effect? 
Cardiology,. 87(4):27 1-275, 1996; incorporated herein by reference). 

While the mechanism by which patients accrue benefit from EECP treatments remains 
unclear, evidence points to the importance of factors other than the obvious mechanical ones that 
prompted early studies of EECP as an external cardiac assist method. It is now thought that the 
benefits of EECP may be related to the recruitment of collateral vessels in the coronary 
circulation, perhaps due to an increase in the synthesis and release of vascular growth factors 
(Soran, A.U., Crawford, L.E., Schneider, V.M., and Feldman, A.M. Enhanced external 
counterpulsation in the management of patients with cardiovascular disease. Clinical Cardiology, 
22(3): 173-178, 1999; incorporated herein by reference). 

Recently, Soran et al have proposed that EECP may act by altering endothelial function 
due to changes in the level of shear stress in the arteries (Soran, A.U., Crawford, L.E., Schneider, 
V.M., and Feldman, A.M. Enhanced external counterpulsation in the management of patients 
with cardiovascular disease. Clinical Cardiology, 22(3): 173-178, 1999; incorporated herein by 
reference). Although EECP has been shown to increase the perfusion through the carotid and 
renal arteries by approximately 20% in one study (Applebaum, R.M., Kasliwal, R., Tunick, P. A., 
Konecky, N., Katz, E., Trehan, N., Kronzon, I. Sequential external counterpulsation increases 
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cerebral and renal blood flow. American Heart Journal, 133(6):61 1-615, 1997; incorporated 
herein by reference), the effects of elevated shear on the synthesis and release of various 
cytokines and growth factors are clearly not restricted to the coronary vascular bed. Shear stress 
is elevated throughout much of the arterial system, especially in the arteries of the lower 
5 extremities by a combination of flow augmentation and reduced arterial cross-sectional areas, 
producing levels of shear stress up to more than four times normal values. These high levels of 
shear occur during both antegrade and retrograde flow (Fig. 5). While potentially damaging to 

=1 the endothelium, these high shear stresses might also provide benefit by stimulating the release 
of shear-induced angiogenic factors from the arterial endothelium of the lower extremities. This 

SfO has not previously been considered and clearly deserves further study. 

f ' : In terms of designing optimal protocols for EECP, it is critically important to understand 

■ 7i the mechanism by which myocardial function is improved. In particular, parameter variations 
\j that optimize the traditional measures of merit (reduction in systolic pressure and increase in 
N diastolic pressure) are not always consonant with the desire to maximize the magnitude and 
15 spatial extent of changes in arterial shear stress in the lower extremity. These issues will need to 
be better understood before a rationale design of the EECP protocol is possible. 

Other Embodiments 

20 Those of ordinary skill in the art will readily appreciate that the foregoing represents 

merely certain preferred embodiments of the invention. Various changes and modifications to 
the procedures and compositions described above can be made without departing from the spirit 
or scope of the present invention, as set forth in the following claims. 
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